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On the dynamics of interfaces in the ferromagnetic XXZ 
chain under weak perturbations 

Bruno Nachtergaele, Wolfgang L Spitzer, and Shannon Starr 

Abstract. We study the time evolution of interfaces of the ferromagnetic 
XXZ chain in a magnetic field. A scaling limit is introduced where the strength 
of the magnetic field tends to zero and the microscopic time to infinity while 
keeping their product constant. The leading term and its first correction are 
determined and further analyzed in more detail for the case of a uniform mag- 
netic field. 



1. Introduction 

The quantum ferromagnetic Heisenberg model has received a great deal of 
interest in recent years. We shall concentrate here on the one-dimensional model. 
The starting point was the discovery of interface ground states in the anisotropic 
XXZ Heisenberg chain by Alcaraz, Salinas and Wreszinski [2] and independently 
by Gottstein and Werner [10]. These interface ground states display two regions 
of opposite magnetization separated by a domain wall whose thickness depends on 
the anisotropy. Explicit expressions for the ground states exist for all values of the 
spin, see Section 2. 

The next important contribution was the proof of completeness of the infinite 
volume ground states (in the sense of local stability) and the proof of a positive 
gap above the ground states for any amount of anisotropy and any value of the 
spin. The spin-i model is SU q (2) symmetric which was used to determine the 
exact value of the gap. This was accomplished by Koma and Nachtergaele [12, 13] 
and Matsui [16]. 

The existence of a spectral gap for general j was proven by Koma, Nachtergaele 
and Starr [14]. Recently, it was proven by Caputo and Martinelli [8] that the gap 
for spin j scales asymptotically like j, however, no precise values are known. This 
verifies a conjecture stated in [14], and comes as a surprise since the spin-j Ising 
chain (which can be considered as the XXZ chain with infinite anisotropy) always 
has a gap equal to 1. 

1991 Mathematics Subject Classification. 82C10;82C24. 

Key words and phrases. XXZ chain, interfaces, dynamics, kink states, Heisenberg ferro- 
magnct. 

This material is based on work supported by the National Science Foundation under Grant 
No. DMS0070774. Shannon Starr is a National Science Foundation Postdoctoral Fellow. 



©2002 by the authors. This article may be reproduced in its entirety for non-commercial purposes. 



2 BRUNO NACHTERGAELE, WOLFGANG L SPITZER, AND SHANNON STARR 

There has been a lot of progress over the last two decades on the dynamics of 
classical spin systems, or lattics gas systems. For instance, for the Ising model the 
Glauber and Kawasaki dynamics have been studied extensively. Here, we only want 
to mention the work by Schonmann and Shlosman [19] on the relaxation in the Ising 
model which is similar to what we want to study in the quantum XXZ chain. A 
quantum model which has been studied in detail is the one-dimensional XY model, 
by Antal, Racz, Rakos, and Schutz [3] at zero temperature and by Ogata [18] 
at non-zero temperature. More recently, Aschbacher and Pillct [5] studied non- 
equilibrium steady states in the same model at non-zero temperature. This is of 
great interest also in the context of entropy production. The XY model (in one 
dimension) is easier to deal with than the XXZ model since it can be mapped to a 
model of free fermions (see [15] and [4]). 

Interfaces and their response to magnetic fields are of great physical and techno- 
logical interest, for instance in storing and reading information on magnetic devices. 
The concept of a magnetic domain was first introduced by Pierre Weiss, already 
in 1906. It was in 1935 in a seminal paper that Landau and Lifshitz presented a 
theoretical treatment of domain walls and their motion when a magnetic field is 
applied. Although the Landau-Lifshitz equation was rigorously derived in the mean- 
field limit by Moser, Prets, and Spitzer [17], its predictions do not universally hold 
true for all magnetic materials. Thus we want to understand the dynamical behav- 
ior of domain walls in a quantum model, in which they occur naturally as ground 
states, namely in the XXZ model. 

The dynamics we consider is generated by the XXZ kink Hamiltonian perturbed 
by a magnetic field. As in many similar cases, we introduce a scaling where the 
strength of the perturbation, say the coupling constant A for a magnetic field, tends 
to and the microscopic time scale to infinity such that their product, which we 
denote by t, is kept fixed. We obtain the following results about this limit: (i) 
For external magnetic fields of bounded spatial support that are sufficiently regular 
functions of time, the leading contribution to the time evolution is identified as 
the reduced dynamics determined by the time averaged field projected onto the 
ground state space of the unperturbed model with an error term of order A 1_<5 , S G 
(0, 1) (Theorem 3.3); (ii) If the spectrum of the unperturbed model has a gap 
between the ground states and the rest of the spectrum, and the rest of the spectrum 
is absolutely continuous, the next-to-leading order term is of order A (Theorem 
3.5); (iii) Analysis of the reduced dynamics for uniform, time-independent magnetic 
fields, reveals a markedly different behavior depending on whether or not the field 
has a non- vanishing component in the z direction. If the z component vanishes, 
the reduced dynamics for the magnetization profile is "ballistic" , while if there is a 
non-vanishing component of the field in the z direction, the magnetization profile 
evolves periodically. 

In Section 2 we define the model and state the assumptions on the magnetic 
field. We derive the leading dynamics and its first correction in Section 3. In 
Section 4 we analyze the leading dynamical behavior in more detail for a uniform 
field and comment on small perturbations thereof. Two appendices provide some 
auxiliary results on the magnetization profile in the kink ground states of the XXZ 
chain and on the spectrum of the Stark-Jacobi operator. 
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2. The model 

2.1. Finite chain Hamiltonian. In this paper we shall only consider the 
spin-i chain. Unless stated otherwise, our results extend to higher values of the 
spin. As usual, let us denote by A > 1 the anisotropy parameter. The kink 
Hamiltonian on the chain [a, b] n Z is defined as 
(1) 



6-1 r 

H[a,b} = -J2 



^VT^A=2(S2-Sf), 



where the spin operators at position x, S X ,S X ,S 3 satisfy the usual commutation 
relations, [S X ,S^] = ie^ 1 SJ. We have added the kink boundary condition, 
— — A~ 2 (Sl — S 3 ) which makes the model SU q {2) invariant and displays the 
full set of kink ground states which we derive next. It is convenient to introduce 
another parameter, < q < 1, such that q + q^ 1 = 2A. One can show that 



6-1 

c 

xx+1 ' 



where P xx+1 1S the orthogonal projection onto the vector g|TI) — lit)- By introduc- 
ing the (j-dcformed lowering operator 



t x +l <S> • • • ® tb 



with t — q 2S ' 3 , we can express all b — a + 2 ground states as (S' _ ) fe |t ■■ ■ T),fe = 

0. 1. . . . , b — a + 1, since [H, S~] = 0. Clearly, these vectors run through all possible 
cigenspaces of the total S 3 component, which obviously commutes with H[a,b]. 
The uniqueness of the ground states follows, for instance, by a Perron-Frobenius 
argument. 

More explicitly, let us denote the eigenvectors of the total S 3, component by 
\(m x )}, where (m x ) £ {±i} b ~ a+1 and Sy\(m x )) — m y \(m x )). Then, as was found 
by Alcaraz, Salinas and Wreszinski [2], these ground states are (up to normaliza- 
tion) 

b 

(2) t=En^ (i_mi)/2 iw), 

(nix) x=a 

where the sum runs over all sets (m x ) such that ^ x m x = m, which is the total 
S 3 eigenvalue. As the ip m are not normalized, we also define \m) = V'm/||V'm||- If 
m = ±(& — a + 1) is maximal/minimal, then we have the all spin-up/down state, 

1. e., the magnetization profile in the z direction equals (m\S x \m) — ±| for all x. 

2.2. Infinite chain kink Hamiltonian. The infinite volume kink Hamilton- 
ian is defined in the standard way via the generator of the Hciscnbcrg dynamics on 
the algebra of quasi-local observables, A. This algebra is the norm closure of the 
algebra of local observables, 

Ai oc = (J A A , A a = (g)Mat c (2). 

AcZ,|A|<oo xeA 
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We prefer to work in the GNS Hilbert space representation of the infinite volume 
kink states. To this end we introduce the incomplete tensor product Hilbert space 



(3) H= |J (g)C 2 ® (g) Q(y) ] , 

A,|A|<oo \xGA j/£Z\A 

where 

IT) ify<o 



im "\\. ir,:>o • 

Let us define the vector 

(4) fi = (g)0(t/), 

which is a vector in H. We also define the (unnormalizcd) GNS vector 

fe>0m<---<£i; t ..<0<iyi<---< 2 / fc j=l 

The set of vectors ip — A£l GNS , where A is a local observable, constitute a dense 
subspace of TL. In this representation, the kink Hamiltonian is now defined by 

(5) HAn GNS = lim [H[a,blA\n GNS . 

[a,b]^Z 

Next we state the assumptions on the magnetic field by which we perturb the 
kink Hamiltonian. 

Assumption 2.1. Let Zxl3 (x,t) ^ B(x,t) e M 3 (ienofe a magnetic field 
which has spatial support in some finite set A C Z, and is absolutely continuous 
with bounded derivative (as a function of time). Let V(t) = 'Yl lxe j j B{x 1 t)S x . We 
further impose a scaling on the strength and the variation of V , and define 

(6) V x {t) = AV(At). 

Clearly H + V x (t) is a selfadjoint operator with domain V{H). The reason 
for our assumption on the support of B is that if we want to consider, for simplic- 
ity, a time independent and translation invariant magnetic field, then, in general, 
e -it{H+v) j g not a strongly continuous group of unitaries on T>{H). The application 
of V on a ground state will no longer be in the Hilbert space TL. This is no con- 
tradiction to the fact that the unperturbed time evolution is strongly continuous 
in the Heisenberg picture, i.e., on the algebra of (quasi-local) observables A, since 
V would no longer be quasi-local. 

We could extend the class in (2.1) to perturbations relatively bounded to H 
(with relative bound less than 1 to ensure selfadjointness according to the Theorem 
by Kato-Rellich) but this does also not allow to consider uniform fields, and we 
therefore stick to bounded perturbations. 

We shall first find the leading dynamics, and then consider the limit as the 
spatial support of B(x, t) tends to infinity. 
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3. Dynamics 

Lemma 3.1 (Dyson series). Let U x (t,to) be the solution to the time dependent 
Schrodinger equation 

ijU x (t,t ) = (H + V x (t))U x (t,t ) 
with V x satisfying (2.1). Then for 4> G V(H) and N e N 
e ltH U x (t,t Q )^ = <f+ H) n / dtV x (t 1 ,t 1 )---V x (t n ,t n )cf> 

0<n<N Jto<t n < — <t-L<t 

(7) +Ejv0, 

where we used the notation alt = dti---dt n , and V(t,t) = e ltH V(t)e~ %tH . The 

error term, Ejv ; is bounded by (A|| V|| ) JV ~ 1 , which converges to as N — > oo 

uniformly for t on compacts; we use \\V\\ = sup{|| V(s)||; s > to}. 

PROOF. We may assume that t = and A = 1. Setting fi(t) = e ttH U(t), then 
for 4> G V(H) we have -^Q(t)4> = —iV(t,t)£l(t)cj>. Integrating over t and iterating 
we obtain 

<„ <Ar ^0<t„<...<ti<t 

+ (-«) Ar / dtV(ti,ti)-"V(*iv_i,tjv_i)f2(tjv-i)^. 

J0<t N -i<...<ti<t 

The error estimate follows then easily. □ 



Figure 1. The first 7 graphs G, comprising Qi, G2 and Qz. 

We discuss next a formula for the iterated t integral allowing us to decide 
which energy terms contribute at what rate. To this end we consider the set Q n of 
directed graphs G with vertex set equal to {n + 1, n, . . . , 1} and bonds satisfying 
the following rules. 1) The direction of a bond is from a higher to a lower value 
of n. 2) There is exactly one bond for every site 1 < k < n, leading to it from 
a site p(k) > k, although many bonds can originate from a certain site. We call 
the bond b(k) = (p(k),k), and define £(b(k)) = p(k) — k the length of this bond, 
and 1(G) — ^£{b(k)) the length of the graph G. 3) If there is a bond between 
i and j, then there is no bond possible between k and I if (k > b(i) and I < 
b(i)) or (k > i and / < i) or (k < b(i), I > i). In other words, there are no 
crossing bonds, and no bonds properly nested. (With these rules, the digraphs G 
are equivalent to the set of all compositions, i.e., all finite ordered sequences of 
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positive integers. To a digraph G = {b(l), . . . , b(n)}, we associate the composition 
[p(l) — l,p(p(l)) — p(l),p(p(p(l))) — p(p(l)), . . .}. Q n is the set of all compositions 
which sum to n.) 

We can show that \Q n \ — 2"~ 1 . It is clear that \Q\\ = 1. If we have a graph 
G (z G n , then we may construct two new graphs in Q n +i- Let us denote the n + 2nd 
vertex by 0, then may be connected to 1 whose graph we denote by G\. The only 
other possibility is to join with p(l) and we call the graph C?2- We define signs for 
the graphs, inductively. Let rx be the unique graph with vertex set {2, 1} and bond 
from 2 (left) to 1 (right), then we set a(rx) = 1. Using the above construction, we 
define the sign by <j(G\) = —a(G), and a(G2) = c(G). Some graphs are shown in 
Figure 1. 

Given a sequence of numbers kj : 1 < j < n, for any graph G € Q n and vertex 
i we define k(i; G) = Ei<j< i+ *(6(i))-i k r 

Lemma 3.2. Let Ej : 1 < j < n + 1 and kj : 1 < j < n be two sequences of 
numbers, then 

n 

^itjiEj+^-Ej+iXkj) _ 




-it(E pW -E 1 +i\k(l;G)) 



PROOF. We can easily check this by induction. The formula is obviously correct 
for n = 1. Assuming that the formula is true for n — 1, we may write the n fold 
integral as 

e -(go (n 



J 1 = i 2 ^o-)-^ + iAfc (i; G ) / 

X y dti ^ e - it i(. E p(2)-E 1 +i\k(2;G)+iXk 1 ) _ e -iti(E 2 -E 1 +i\ki 

The vertex set of the graphs G £ Q n -i is here {n+1, . . . , 2}. Now, in the first term of 
the integrand we have that p(2; G) = p(l; H), k(l; H) = k(2;G) + k 1 ,a(H) = a(G), 
where H is the graph on {n + 1, . . . , 1} which agrees on n + 1, . . . , 2 with G plus 
a bond connecting 1 with the same site as 2. In the second term we have a bond 
between 1 and 2, and by our rules hence a change of the sign. We can combine this 
into a sum over all graphs H e Q n , which verifies the formula. □ 

In the following theorem we identify the leading term in the limit when A — > 
0, and t —> oo, such that r = Xt is constant. Here, r can be interpreted as a 
macroscopic time scale. We shall use the spectral decomposition of the Hamiltonian 



(9) H= EdP(E) 

Jo 

such that for any <j> £ H, we have 

r-OO 

(j)= dP{E)(j> . 
Jo 
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Theorem 3.3 (Leading dynamics). Let <j> S Ti., and let V satisfy Assumption 
(2.1). Then, for all 5, 1 > 5 > 0, there exists a Ao depending on \\V\\,t, \\(j)\\ and S 
such that for < A < A 



(10) 



POO 

^U\X-^-j T (e-iP(E)5\]dtv{t)P{E)^ dp{E)(j) 



< A 



1-5 



T means time ordering, i.e., T(V(s)V(t)) = V(s)V(t) if s > t, and zero otherwise. 

PROOF. We may assume that \\<fi\\ = 1. Let N be an integer which we choose 
later. We may write (t = A _1 r) 

JtH^^-i^ - J°°T (e- lP ^ lo dtv(t)P(E)^ dp{E)(t) 

< \E N \ + |Ejv - J? N \ + \E' N \ . 

Here, E^v is the error term from equation (7), and the remainder term of the 

power series for the exponential function in T (e~ tP ( E ^ ■A) dt y (*) p (- E )^ <fi to the order 

N — 1, whereas E^v and E'^, refers to their finite sums. 

First we take care of the error terms E^ and ~E' N by choosing N = N(X;\\V\\,t) 
such that 

1 (r| 



N 



1 \ 1-5 



\\V\\ N\ <3 A • 
This can be accomplished with TV = — Cln A and the constant C depending on ||V|| 
and t. 

Second, we investigate the limit (A — > 0, t — > oo, At — > r) for the n-th term 
in the finite Dyson series. By inserting the spectral resolution for H (9), the n-th 
term in e ltH U x (t)(j> equals 



(ii) h 



A)" / dt f dE n f[ e -itAE J+1 - Ej ) T I TT p(Ej)V(\tj) I dE n+1 4> . 



In order to perform the t integration we write V as its inverse Laplace transform. 
Let ci < C2 < • • • < c„ < 0, then we obtain 

rci+ioo />c n +ioo 

dk 



HA)"— — / •••/ dk dE n 

(2m) n J Cl _ loo J Cn -ioo Jr™ +1 

« n In 

x dt TT e - it *( E *+ 1 - E ' +iXk ri T ff PtEJVfc) | dP(E n+1 )<f>. 



By using formula (8) we can do the t integration, and obtain 

i /-ci+zoo />c n -\-ioo f- f- n 

— V / ■■■/ dk dE dtT] e- lt ^ E ^- E ^' 

(2m) n J Cl - loo J Cn - toa J R ^+! y [0 ,t]« fJl 



e °& n 



ceg„ 



-it(E p{1) -E!+iXk(l;G)) _ 
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Here, it is convenient to have q's such that there are no zeros in the denominator. 
We observe here the crucial fact that in order to find a term of the order 1 , all Ej 
have to be equal. We may now go back to (11), change coordinates and conclude 
that in this limit we get 



HA) 



"/ / rftT [ TT P(E)V(Xtj) I dP(E)(j> 
Jo J[o,t]» \fJi J 

= {-i) n J r(J dtP{E)V{t)\ dP(E)<p. 

Finally, we have to bound the sum of all other terms in the finite Dyson series 
which individually tend to as A — > 0, and are contained in — ^jvl- As there 
are 2 n_1 terms in formula (8), we immediately obtain the bound 

A^ 1 (2r||y||)" 
\L, N ~ L N \ < - 2 , 



2 ^ n\ 

n—l 

which in turn is bounded by |A 1_(5 if we choose A < A with A such that e 2r " y " < 

2 \~ s n 

3 A • u 

If H(f> = is a kink ground state of H , then the leading term equals 
(12) T^/^MWojj^ 

This expression will be further analyzed for a specific choice of V in Section 4. 

Note that so far we did not use any specific spectral properties of the Hamil- 
tonian, and thus Theorem 3.3 is valid, for instance, for any value of the spin. 

We know that <j(H) C {0} U [1 — A -1 , oo), using the theorem of [12] that there 
is a gap of 1 — A -1 above the infinite volume ground states. If we decompose 
the kink Hamiltonian on Z as the sum of a kink plus a droplet Hamiltonian, i.e., 
= H^Soo y ] + H~[y~oo) w ith y G Z, then we get approximate eigenvectors by 
tensoring a kink ground on (— oo,y] (and with interface far to the left of y) with 
eigenvectors of HT^y In the case y = — oo the spectrum of the latter has been 
shown by Babbitt and Gutkin [6] to be {0} U [1 — A _1 ,oo). Their analysis follows 
the earlier work by Babbitt and Thomas [7] for A = 1 proving the completeness 
of the Bethe ansatz eigenfunctions in the infinite volume limit and thereby also 
showing that the spectrum is absolutely continuous. 

This way, one obtains a(H) — {0} U [1 — A -1 , oo). But this construction does 
not give the nature of the spectrum in the interval in [1 — A _1 ,oo). We plan to 
address this issue in the future; for now we state the following hypothesis. 

Hypothesis 3.4. The spectrum of the spin-\ kink Hamiltonian H with A > 1 
has no other point spectrum than 0, and the rest is absolutely continuous, equal to 
[l-A-\co). 

By assuming this information we can identify the first correction to the leading 
term given in Theorem 3.3 for ground states, and we show that it is 0(A). For 
simplicity, we state the result only for time independent fields and leave the general 
case to the reader. The main argument in the proof below, however, is carried out 
for a time dependent field. 
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Theorem 3.5 (First order correction). Let <f> be a (normalized) ground state 
of the spin- j kink Hamiltonian, V a perturbation as in Assumption (2.1) and H 
satisfying Hypothesis (3.4-). Then, 

giA-^Hg-iA-MJT+AV^ _ e -irP(0) V(j) _ x T ^p {E)Ve ^rP(E)V ^ 

Ji-a-i E 

(13) - P(0)V (e-iT[P(o)v+AH-i(i-p m v] _ e -irP(0)v^ J = o(A) _ 

Proof. We analyze the n-th term in the finite Dyson series and treat the error 
terms as before by choosing first N = N(\, ||^||,r). The spectral resolution of H 
has now the form 1 = P(0) + f™ A -i dP(E). We have seen above that to leading 
order all energies Ei have to be equal to 0. According to (8), terms of the order 
A come from sequences (E m — E{) : {m < 1} € G such that all but one pair is 
non-zero. We represent such sequences graphically by polygons with n + 1 vertices. 
We draw a horizontal line between the vertices k and k + 1 whenever Ek = Pfc+i, 
an upward line if Ek > Ek+i = 0, and a downward line if = Ek < E^+i- 

By our rules for the graphs G e Q n , the only possible polygons contributing to 
the order 0(A) are of the form displayed in Figure 2. 




Figure 2. Sequences of energies satisfying the condition that 
their corresponding differences is non-zero for exactly one pair. 
The bump in polygon P 2 can be at any site such that the starting 
and ending point of the polygon lies at the horizontal line. 



Using Hypothesis 3.4 we shall show that, in fact, only Pi and P2 contribute to 
0(A) while P 3 and P4 are of lower order. We shall now consider the contributions 
from Pi, P 2 , P3 and P4 to formula (11). Let us start with polygon Pi, i.e., we have 
a sequence = E n+1 < E n = . . . = E\ = E. In 

rtn-l 

/ dt n e lt - E V{\t n )4> = 
Jo 

-L (e^-^ViXt^) - V(0)) <f>-± J^" dt n e u ~ E V{\t n )<j> 

we can use that, by our assumption on the spectrum of H and the Riemann- 
Lebesgue Lemma, the leading term is -^V(0), and the rest is o(l), where o(l) 
depends on the inverse gap A/(A — 1) and sup{||V r (s)|| : s e [0, r]}. Then, up to 
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o(l) we have 

r°° 1 r rc— 1 

a / dE -Hxy 1 - 1 / dt n p(^)y(Ai,)P(£;)y(o)^ 

dE-H)"" 1 !^ dnP(£;)y(Ti)J P(^)F(O)0, 

where we have used a change of variables as in the proof of Theorem 3.3. This gives 
the first correction term in (13). 

Next, we consider the polygon P2 with all energies equal to except for E m = 
E > with 2 < m < n. Then, the contribution of P2 to the n-th term in the Dyson 
series is equal to 



{-i\) m+1 / dE / dt m e~ u ™- lE+ltmE 

Jl-A- 1 Jo<t m <---<t!<t 



m— 1 

n 



A(Ai m ) 



with 



^ n 

A(\t m )=P{E)V{\t m ){-iX) n - m - 1 dt n _ m Y[ P{0)V{Xtj)<l>. 

J0<t n <-<t m j=m+1 



Note that A(Q) — 0. In order to simplify matters, we let V be time independent 
and perform the time integration. Hence, we get (to leading order) 

( — iAt) n_1 f°° 1 / \ro-l / \r 



(— i\t) n ~ l f 1 / \m-l / 

= a t ttt(P(0)^) ff _1 (l-.P(0))V (P(O)V 



(n- 1)! 

Now we can easily sum up the contributions of all n — 2 polygons of the form P2, and 
obtain the n-th term in the power series of the exponential in the second correction 
term in (13) up to errors o(A). 

Finally, we prove that the other two kinds of polygons contribute of lower order. 
For the polygon P 3 , (11) equals 

ft 



(-iA) / dt 1 e ltlE P(E)B(Xt 1 ) 
Jo 



with 



B(Xti) = (-iA) n_1 / dtV(Ati) JJP(O)V(Atj)0- 

J 0<t n < - <tx j =2 

B(0) = 0, and by integration by parts we get 

ft \ 
(-iA) / dh e ttlE P(E)B(Xt 1 ) = --e ttE P(E)B(Xt) + o(A) 
Jo E 

by our assumption on H . 

For a sequence of the form P4, with E = E\ = . . . = E m > E m+ \ = 
E n +i = 0. Then, the expression in formula (11) equals 



(-iA)' 



r. m 

■ dte i *- B TTP( J B)y(At j )C(Af ro ), 

Jo<t m <-<t - =1 
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with 

c(\t m ) = {-iX) n ~ m / dt TT p(o)v(xtj)4>. 

Jo<t n <-<t m j=m+1 

Integration by parts in the first integral shows that this term is o(A); in fact, if we 
assumed V to be m times diffcrcntiablc, then it would be o(A m ). □ 

Notice also that the first correction term of the order A contains the information 
on the whole spectrum while the leading term only depends on the highly reduced 
space of ground states, which we analyze in the next Section. 



4. Leading dynamics 

In this Section we consider the leading dynamics generated by P(0)VP(0). We 
start with a uniform magnetic field V\ for which we determine an explicit expression 
as A — > oo. 

4.1. Uniform time-independent magnetic field. Let B(x,t) = B be a 
constant vector, and thus V A = X^eA B ■ S x . We shall consider K\ — P(0)V\P(0) 
as a Jacobi operator (i.e., tridiagonal) on £ 2 (Z). Since H from (1) commutes with 
global rotations in the xy plane, we choose a new basis \n, U) = ?7 _1 |n) of P(0)H, 
where U satisfies UVaU' 1 = Y,xeA V B 'f + B^S* + B 3 S%. In other words, U = 
l\ x e I0S *, where {B U B 2 ) = ( B{ + B% cos 9, y/Bf + B% sin 9) . 

When taking the infinite volume limit we have to renormalize the total S* 3 
component, and define 

(14) S 3 = lim ]T S 3 X . 

\x\<n 

In the new basis, we have 

(m,U\K A \n,U) 

= (m\P(0)UV A U- l P(0)\n) 

= + B*(m\[S+ + 5-] |n) + B 3 {m\S%\n) 

xeA 

= ^Bl + Bl o A (n) (S(m, n - 1) + 5(m, n + 1)) + B 3 b A (n)S(n, m) , 

where the index n G Z denotes the eigenvalue of the renormalized total S 3 compo- 
nent. In Appendix A we calculate a = limA^zaA^), see equation (22). Notice 
that the limit is independent of n by translation invariance. The function 6a (n) 
satisfies |6a(ti)| < |A|, and limA^z b\(n) — n for all n. 

For finite A, K\ is a trace class operator, and the spectrum is pure point, 
accumulating to zero. We now discuss the case when A = Z. 

Let A be the lattice Laplacian with unit off-diagonal entries and zero diagonal 
terms on £ 2 (Z), and let W(n) = n be the Stark potential. Then we define the 
Stark- Jacobi operator, 

(15) K Q = aA + ~/W 
with a = \/B\ + B\a and 7 = B 3 . 
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The operator Kq is selfadjoint with dense domain T>(Kq) = {/ e £ 2 (Z) : 
(nf(n)) n e £ 2 (Z)}. This follows from the fact that W is selfadjoint on T>(Kq), and 
that A is bounded and selfadjoint. 

As for the continuous case, one can find the spectrum explicitly but it turns 
out to have rather different properties. Only in the case 7 = we shall have pure 
absolutely continuous spectrum which is equal to [—2a, 2a]. As soon as 7 ^ 0, the 
spectrum is pure point. The latter is analyzed in the following lemma, and the 
generalization to higher dimensions is given in Appendix B. 

Lemma 4.1. Let Kq = aA + -fW be the Stark- Jacobi operator (15), and 7 7^ 0. 
Then, spec(Ko) = 7Z is pure point spectrum with simple eigenvectors 4> m (n) — 
J rn - n (2a/"f); J n denotes the (ordinary) Bess el function of order n. 

PROOF. By Fourier transforming the eigenvalue equation a<f>{n — 1) + a<p(n + 
1) + -fncf>(n) = \<p{n), we obtain 

[ 2a cos k + 17-77- ) 4>{k) — X(j)(k), 
\ dk J 

where k e [— K,ir]. This is solved by cf>(k) — e *( 2 <*sm/c-Afc)/7. The periodicity 
condition cf){k + 2n) = (j)(k) shows that A G 7Z. Note, that X^mez = 2n8(k) 

which proves that these (mutually orthogonal) eigenfunctions are complete. In 
configuration space we thus have 

Mn) = ^/ cos(^sinfc + /c(n-m))dfc = J„- m {^) ■ 

□ 

Lemma 4.2. The kernel for the reduced dynamics generated by P(0)VP(0), and 
hence of K n , has the following form: 

(x,U\e~ up ^ p ^\n,U) 
(16) 

Proof. 

e- tKot S n (x) 



The last step is content of Graf's addition theorem, see [1], formula (9.1.79). □ 

For future ease of notation, we define w = — sm ("jf); an d X = ■ 
We now inspect closer the reduced time evolution of the magnetization profile 
of a kink ground state localized at 0. Defining 

m(Cl, x, t) = (0, U\e up ^ VP ^ (Cl ■ S x )e- ltp ^ vp ^\0, U) 



(x) 



= J„_ x (^ sin (^))ex P U( 



"yt — 71 -ft + 71 



„)■ 



^ mt ^ m (x)^ m \Sn) 

E -i^mt j (2a\ j (2a\ 

E -i-ynt i-ymt j i la \ j /2a \ 

e c •Jm — n+x\ j )^m\ ^ ) 



e -ijnt e i(w-jt)(x-n)/2 



J n - X (>sin(f; 
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for any vector ft € S 2 , we calculate 
m(0, x, t) 

= Y, < ' U\e up ^ VP ^\m, U) (m, U\fl ■ S 3 \m' , U) (m' , u\e' up ^ VP ^ |0, U) 

m.m'GZ 

= ]T e ix(ro '- m) ^H^HH«-e(fi)-5 x |m'), 

m.m' EZ 

where is rotation about the z axis by an angle —6. (Note this is the opposite 
rotation as that used to go from P(0)VP(0) to K .) Taking = £3, and defining 

m 3 (x,t) — m(es,x,t), we have 

- 3 (s,t) = Yl e^ m '- m h m (w)J m ,(w)(m\S 3 x \m'} 



m 



Yj 2 m (^sm(%))( m \S 3 \m). 



Similarly, defining m 1 (x,t) = m(ei,x,t), and observing that R_e{e\) — ae\ — be 2l 
where a = cos 6 = B\ / y 1 B\ + B\ and b = sin = B 2 / y 1 B\ + Sf, we have 



^ >->ir n ^ 



|m'). 



We now observe a phase condition on |m) that guarantees (mlS 1 |m') is real. This 
means that 



m 1 (x, t) 

= e^ m '- m \j m (w)J„Aw) 



e lx (a-ib) e~ ix (a + ib) 



a — ib 



, a + ib , 



, {m\S x \m!) + -^—{m'\S x \m) 



Y J m {w)J m +i(w)(m\S x \m+ 1) 



nit; 



2 2 
= [acos(x) + 6sin(x)] Y Jm{w)J m+1 {w){m\S x \m+l) 

mGZ 

= [ as in(f) + 6cos(f)]^ J ro (4fsin(f))J m+1 (4fsin(f))(m|5-|m+l). 

m£Z 

Notice that for 7^0 the motion is periodic whose period is solely 

determined by the field in the z direction. The profile m 3 (x,t) is given by an 
absolutely convergent series. Therefore, pointwise convergence implies convergence 
on compacts. As the function is periodic in t and thus the image is a compact 
interval, the profile stays exponentially localized for all times. 

However, as 7 — > 

(17) m 3 =0 (x,t) = Y Jl(2at){m\Sl\m), 

mGZ 

(18) m} =0 (s,t) - bY J m (2at)J m+1 (2at) (m\S+\m - 1) . 

mGZ 
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Next, we analyze m 3 (x,t) in more detail. A very similar situation has been 
studied by Antal et al [3] in the context of the time evolution of the XX model. 
We follow their analysis. Although we have just derived explicit formulas for the 
profiles, which can be analyzed numerically, we prefer to have simpler expressions. 
As we shall show, for large x and t, m 3 (x, t) will be a function of the velocity v = j , 
only. To this end we define the discrete derivative of m 3 (x,t) as 

<j>' x (v) := t\m* (x + l,t)-m* (x,t) 

J f=v 

and study the limit lim^oo (j)' x {v). Note that (0\S 3 _ m _ 1 - S^_ m |0) > for all 
x and m because the profile of the kinks is a decreasing function. If we define 
p(m) — (0|<S^_ 1 — S'mlO), then since it is a telescoping sum 

oo 

]T p(m) - lim (0|^|0) - lim (0|S» = 1 , 

*• — ' m— >oc m— *— oo 

m— — oo 

i.e., p(m) is a probability distribution. All moments of p are finite and the first 
moment of p is |. In fact p(m) <~ e~ c l m for some positive c depending on A > 1 
(see Lemma A. 1.1). Thus, 

C(«) = -Ef J ™+*(irMm). 

mGZ 

We used the fact that J_„(x) 2 = J n (—x) 2 = J n (x) 2 . One may easily bound 
v *^r-m(ir) uniformly in x and m, which our asymptotics will show (see below). 
Hence, we may interchange the limits x — > oo and the sum, to calculate 

lim = ]T p(m) lim fJ^ +x (^) . 

We distinguish between the two cases v > 2a, and v < 2a. The asymptotic 
behavior will be very different in these two cases. Let us start with v > 2a. Then 
by (9.3.2) in [1], asymptotically for large x 

2 2£>J* (2*°) = 2(m+x)a j2 / / + ) 2ga \ 

v m+x\ v I m+x x m+x \ \ ' (m+xjv J 

1 _2as(i_J3^ZI) 

r^j P v V V 4a 2 1 

Since p(m) is a probability measure, which decays exponentially, the bound above 
shows that 4>' x {v) is exponentially small, and lim^^oo 4>' x {v) = 0, which implies that 
(p(v) := lim^^oo 4> x {v) = —\. This behavior is not surprising for it just means a 
finite speed of propagation, which has to be less than 2a. 

However for < v < 2a, the behavior of the Bessel functions changes, cf (9.3.3), 

[1]: 

2 -^J 2 m+ x( 2 -T) ~ , 2 2 cos 2 (z^-l-zaxccos (£)-£) . 
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This implies that for large x 

<t>' x {v) ~ jL== cos 2 [x^-l-x arccos (£) - f ) . 



an^l 

For v < 2a we recover the function 4>{v) by integration, (we switch the limit and 
integral by dominated convergence, since (f> x {v) is uniformly bounded in x and v 
iorve [0,|-e]), i.e., 



(j)(v) := lim (f) x (v) 



[ v 2 

/ dy lim 4>'(y) = arcsin(^-). 

JO ' x ^°° TT Za 



The oscillating cosine contributes an average factor | to the integral. Using sym- 
metry, and collecting terms, we have shown that 

+\ for v < -2a 



(19) lim_ m„ (=0 (x,t) = (j>(v) = < -f arcsin(^) for -2a < v < 2a 



i for v > 2a 



The importance behind a further analysis of the profile perpendicular to the 
z direction is that it allows us to decide whether the state is rotated by the per- 
turbation into the xy plane or whether there is ballistic diffusion of the interface. 
Comparing equations (17) and (18), we see that the main difference is a replace- 
ment of the probability measure p(m) = (0\Sf n _ 1 — 5^J0) by the signed measure 
p(m) = (l\S+ — S^ l _ 1 \0). One can easily determine that |p(m)| < Cie~ c l m (see 
Lemma A. 1.2). However, now p{m) oscillates, and again because of a telescoping 
sum, 

oo 

]T p(m)= lim(l|S+|0>- lim (1|S+|0> = - = 0. 

m— — oo 

Since we have determined that for any fixed finite m, lmx^oo ^J^ n +x(^r) 1S a 
fixed function of v/2a independent of m, and in particular this implies the same 
limit for lim^oo ^J m+x (^)J m+x+1 (^-), tlie oscillation of p implies that 



lim ip' x (v) := t ml =0 (x + l,t) - ml =0 (x,t) 

x^oo.x—vt L J y~ v 

= lim iVp(m) lim ^J m+x (^)J m+x+1 (^) 

mEZ 

= 0. 

This does not say that the profile in the xy plane is vanishing. Quite on the contrary, 
it means that it stays exponentially concentrated around the kink at 0, and does 
not follow the ballistic motion of the third component. 

Let us shortly comment on time dependent uniform fields with time dependent 
Stark- Jacobi operator K{t) = a(t)A + j(t)W. It is easy to see that if Bs(t) = 
then we just have to replace at by J* a(s) ds. E.g., 

m 3 (x,t) = ^ Jl(2f*a(s)ds) {m\Sl\m) . 



We have not been able to derive explicit expressions in the general case, but one 
can show, for instance, by repeated use of Graf's addition theorem, that there is a 
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function F such that 

m 3 (x,t) = Jl(F(a(t)Mt)){m\SZW- 

In these time dependent cases, a scaling limit as before does not exist. 

4.2. General infinitely extended fields. Since a thorough discussion would 
be out of the scope of this paper, we only want to convince the reader that the 
leading dynamics of the "many-body" Hamiltonian leads to a well studied one- 
body problem which has seen great progress in recent years. 

Let us start with the simplest case when the magnetic field B(n) is asymptoti- 
cally uniform with non-zero third component. I.e., we assume that there is a vector 
B and some p > 1 such that (B 3 (n) - B 3 )n e t?(Z), and (B h2 (n) - B li2 ) G £ P (Z). 
Let K denote the Jacobi operator corresponding to this vector field. Then K is a 
compact perturbation of the Stark- Jacobi operator K in (15). By Weyl's Theorem, 
K has pure point spectrum and the motion is quasi-periodic and the interface stays 
localized at 0. 

The case of a perturbation of a uniform field with _B 3 = is much more com- 
plicated. If we perturb by a field in the z direction, which corresponds to adding a 
potential to the (discrete) Laplace operator, then we are in the widely studied case 
of a (discrete) Schrodinger operator. 

The effect of a perturbation of B in the xy plane, however, is non-local. I.e., if 
we change B\ at into Bi{0), then the off-diagonal matrix element computed from 

(n\ ]T B^S+ln - 1) = BM ]T S+\n - 1) + (B^O) - Bi)(n|S+|n - I) 

is effected for all n, although exponentially localized at the impurity at 0. By a 
rotation, one can map this Jacobi operator into a Schrodinger operator. Any non- 
uniform perturbation of the uniform case by a bounded perturbation will modify 
the spectrum of the reduced Jacobi operator, as was recently proven by Killip and 
Simon [11]. 

There are many more interesting fields for which we would like to understand 
the time evolution of interfaces. A particularly simple case is a sharply localized 
field at a single site y, i.e., V = B ■ S y . Here, one should rather approximate the 
time evolution by the projection of H + B ■ S y onto its low energy spectrum which 
seems to be separated by a gap from the rest of the spectrum, and which is similar 
to the spectrum of the unperturbed kink Hamiltonian. So far only the ground 
state and the gap above the ground state of this Hamiltonian have been studied by 
Contucci, Nachtergaele and Spitzer [9]. We plan to pursue this in the future. 

Appendix A. Magnetization profiles 

The exact formula for the spin-^ magnetization profile of kink states in the 
z direction was derived in [20], pp 55-58. Note that by rotational symmetry, 
(n\Sl\n) — (n\Sx\n) = 0. We introduce the function 

(20) f(z) = ^{-l) k z k q k{k - 1 l 

fc>0 

Lemma A.l. For the spin-^ infinite chain kink Hamiltonian we have the fol- 
lowing profiles of the normalized kink state, |0) (which has total z component 0): 
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(1) lfx>0, then 

1 OO 

(21) (0|S 3 |0) = - - + q 2x J2(-l) k q k{k+2x+1) , 

k=0 

while for x < we have (0|S 3 |0) = -(0\Sf_ x \0). Further, p(m) = 
(0|S'^ l _ 1 - S? n \0) < e _c l" l l with c depending on q. 

(2) (n\S \n - 1) = gM./VM+ 2 ) ; and p(m) = - S+^O) satisfies 
\p(m)\ < Cie _c l" l l with c depending on q. 

(3) The number a in the Stark- Jacobi operator (15) is equal to 

(22) a = <n| £ Sl\n - 1> = \ £(-l)V<^>I±* 

sez fe>o 

Proof. Here, we only show the key steps for Part 2 and 3; the details can be 
found in [20]. The exponential localization of the measures p(m) and p~(m) can be 
deduced from the explicit formulas but it can also be done much easier by using 
the following estimate 

(23) (Mh+*&(*)Sl]\^)<CqW 

for some C = C(q) independent of x. Let S1 H = &) x<0 (J) ®® x >o (?)• We define 
the so-called grand-canonical states 

oo 

(24) i,(z)= J] (l + z- 1 q- x S-)H(l + zq x S+)n+-. 

x— — oo X—l 

The reason for this definition is that 

(25) = 

n6Z 

It is clear that in the grand canonical ground state, for real z = q~^, 

(lP{ Z )\S^(z)) _ 1 g (*-M)/2 _ g (M-*)/2 
(V'WIV'O)) 2 + q ^-x 

which proves exponential localization of the interface in the grand canonical ground 
state. E.g., setting fi = 0, one has for x > 0, 



1 < 1 q x ' 2 - q- x l 2 _ 11 v/TT^ - 1 + 



2~2 ^ + g -* 2 2 ^1 + ^ 

(26) <_! + ![!+ 1,2, + 

and a similar inequality holds for x < 0. Setting fi = 0, (i.e. z = 1) we see that 
(^(l)|(i ± S 3 )|^(l)) = X>*l(2 ± ^ 3 )|^„) > (Vo|(| ± ^ 3 )|^o) , 

using the fact that (| ± S x ) is a nonnegative operator. Hence, 

(M^±^M H(z)\\ 2 Mi)|[j±sgMi)) 
(V>olV>o> - \\M 2 <^(i)|V(i)> 
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Since (||V'( 2; )l| 2 /IIV'o|| 2 ) < oo, and using equation (26) we derive equation (23). This 
proves Part 1. 

Next, using the g-binomial product theorem we get that 

— 1 oo 

(^( Z )\s+mz)) = z- 1 n (i+2-v 2x )n( i+ *v x ) 

X— — oo x — 1 

( 00 \ / 00 „n(n+l) \ 

_ x ~ |2n g |n|(|n| + l) /(g 2(|n| + l) ) 

~ Z l z l (a 2 - a 2 ) 

n=-oo l<7 ,<7 Joe 

Comparing powers of z in (25) we have 
(n\S+\n-l) = 

Un-l\\Un\\ 

= gH/(<7 2 H+ 2 ). 
We have used that ||i/>„|| 2 = q n(n+1) / {q 2 ]q 2 )oo- Finally, 
5>|S+|n-l) = ^(x\S+\x-l) 

= qM ^(-l) fe 9 2(|;r| + 1) ' £+fe(fe " 1) 

x& k>0 

i i „l+2fc 
- Vf 1)k n k(k+l) 1 + g 

^ * 1 _ „l+2fc • 

fe>0 y 



□ 



Appendix B. Stark- Jacobi Operator on Z 

Lemma B.l. Let Kf(n) = (aA + j-n)f(n) be selfadjointly defined onV{K) = 
{/ € ^(Z d ) : (7 • n)f(n) G ^ 2 (Z d )}. // aZZ 7j are non-zero and there is some 
70 smc/i t/iai 7j = aj7o aj G Z /or a// j = 1, . . . , d, i.e. if (71, . . . , 7<j) are 
commensurable, then 

spec(K)= 1171,2 — Z, 

lcm(ai, . . . ,a d )7o 

w/iere lcm(ai, . . . , a d ) is the least common multiple. If all gj are nonzero, but they 
are not commensurable, i.e., there exists a pair ^j, -fk such that^j/^k is irrational, 
then spec(K) = R and the spectrum is dense pure point spectrum. In this case there 
is an eigenvalue at every point of the lattice 

J2^fz = {'£a j \\ 1 \\ 2 / lj :a 1 ,...,a d eZ}. 
3=1 ' 3 3=1 

However, if there are 1 < I < d non-zero components 0/7, say the first I, then 

spec(K) = cl([-2a(d - l),2a(d - l)\ + ^(M^Z)) , 

j=i lj 
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where cl(.) is norm closure in R, and there is essential spectrum even when all jj 
are commensurable. 

Proof. Let us assume first that all components are non-zero. In Fourier space 
we have the eigenvalue equation 

2a £ cos ( kj ) + j|| 7 || £ ^ A - \j m = . 

Let R = R(j) be an orthogonal matrix such that Y^j=i ^y'lpylf = ^(*'^)- (^- n 
particular this means Rij — p^-) Setting k — Rk, we get 

d i ( d d \ 

9fc i nil \ 1=i ' J 

The solution is thus 

cf)(k) = Ce~ ifw( A ^ 1_2Q ^j d =i(- Rl ^ ^"E^i^j^)) _ 
Using that Rij — -pL^ we have 

(27) 0(fc) =(7 e -lBTr( A S 3 d =i 1 ^ r ^-2aE 3 d = i J ^ !1 sin(/ S3 )) ^ 

Finally, the periodicity condition, 0(fcj + 2ir) = <j>(kj),j = 1, . . . , d requires that 

\ e Jh£ z . 

7j 

Now suppose that the first ^ components of 7 are non-zero and the rest are zero. 
According to the decomposition £ 2 (Z d ) = £ 2 (Z') £ 2 (Z d ~ l ), we write if as the sum 
of Stark- Jacobi operator acting on £ 2 (Z l ) and aA^ acting on £ 2 (Z d ~ l ). Hence 
the result follows. □ 
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